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Compact stars made of quark matter rather than confined hadronic matter, are expected to form a 
color superconductor. This superconductor ought to be threaded with rotational vortex lines, within 
which the star's interior magnetic field is at least partially confined. The vortices (and thus magnetic 
flux) would be expelled from the star during stellar spin-down, leading to magnetic reconnection at 
the surface of the star and the prolific production of thermal energy. In this paper, we show that 
this energy release can re-heat quark stars to exceptionally high temperatures, such as observed for 
Soft Gamma Repeaters (SGRs), Anomalous X-Ray pulsars (AXPs), and X-ray dim isolated neutron 
stars (XDINs). Moreover, our numerical investigations of the temperature evolution, spin-down rate, 
and magnetic field behavior of such superconducting quark stars suggest that SGRs, AXPs, and 
XDINs may be linked ancestrally. Finally, we discuss the possibility of a time delay before the star 
enters the color superconducting phase, which can be used to estimate the density at which quarks 
deconfine. From observations, we find this density to be of the order of five times that of nuclear 
saturation. 



I. INTRODUCTION 

The interiors of compact stars provide a naturally- 
occurring environment that can be used for studying the 
properties of ultra-dense baryonic matter. A common 
means of probing this environment is through direct ob- 
servations of thermal emission from the surface of com- 
pact stars. By comparing these observations with the- 
oretical models one can retrieve key information about 
the physical processes occurring in matter compressed to 
ultra- high nuclear densities [3-|g|. 

Most attempts to model thermal emission make use 
of what is called the minimal cooling scenario, which in- 
volves cooling through the minimum set of particle pro- 
cesses that are necessary to explain the thermal evolution 
of the majority of compact stars. Appropriately, it is 
used as a benchmark for observations of cooling neutron 
stars. However, there are a number of compact stars 
(see Table HJ possessing thermal emissions significantly 
out of agreement with the minimal cooling scenario, in- 
dicating that other processes may be occurring and need 
consideration. Some specific classes of compact stars that 
disagree with minimal cooling, and are already distinct 
for some of their other features, are Soft Gamma-ray Re- 
peaters (SGRs), Anomalous X-ray Pulsars (AXPs), and 
X-ray Dim Isolated Neutron stars (XDINs). It is gener- 
ally accepted that SGRs and AXPs are the same type of 
objects, and it has been speculated before that XDINs 
are also related 0j- 

Observations indicate that SGRs and AXPs are very 
hot objects. Heating by magnetic field decay in the crust 
has been suggested [8j as a possible explanation and, 
at certain times during the star's evolution, it has been 
shown to possibly be the dominant source of heating de- 



pending on the values assumed for magnetic field decay 
timescales [9J. However, even with the most liberal de- 
cay timescale parameters, crustal magnetic field decay 
does not produce enough heat to account for SGR and 
AXP observations over the span of their lifetimes 
In fact, there is no (micro-)physical model that can ex- 
plain the temperature evolution of the objects in Table 
|TJ Phenomenological studies have been performed [TTI j . 
which parametrize the magnetic field decay to fit obser- 
vations, and conclude that "some efficient mechanism of 
magnetic flux expulsion from the star's core is required" . 
Other studies proceed by introducing an artificial heat 
source in an internal layer [l2| . 

In this work, we expand on the idea that magnetic field 
decay from compact stars causes a re-heating of such 
objects. As stated above, the net effect of this mech- 
anism on the temperature of ordinary neutron stars is 
much too weak to accommodate the temperatures ob- 
served for SGRs, AXPs, and XDINs. The situation 
changes dramatically if one assumes that these objects 
are made of superconducting quark matter rather than 
confined hadronic matter. As demonstrated in our pre- 
liminary work fl3| and expanded upon herein, magnetic 
flux expulsion from the cores of such stars provides a very 
efficient and robust mechanism that can re-heat com- 
pact stars to the temperature regime observed for AXPs, 
SGRs, and XDINs. 



II. THE VORTEX EXPULSION MECHANISM 

The feature that quark matter forms a (single species) 
color superconductor is critical in our study. The conden- 
sation pattern that is considered here is the color-flavor- 
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TABLE I: Compact stars with reported thermal emissions sig- 
nificantly higher than can be predicted by the minimal cooling 
scenario, specifically SGRs, AXPs, and XDINs [9|. The listed 
ages are the spin-down ages modified by the vortex expulsion 
process (P/3P), instead of the usually assumed spin-down age 
(P/2P). 
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1 1E2259+586 and 4U0141+615 are not considered here. As ar- 
gued in [T3 | these harbor an accreting ring that explains their ex- 
treme luminosities. 



locked (CFL) phase fl5l - fl8l ] where quarks of all colors and 
flavors pair together to form Cooper pairs. Henceforth 
we refer to such stars as color-flavor-locked quark stars 
(CFLQS) [HI. We note that, while in this paper we have 
focused on the spin-zero CFL phase, other color super- 
conducting phases may be more applicable, for instance 
the spin-1 color superconducting phases Because 
of stellar rotation, CFLQSs develop rotationally-induced 
vortices [2(|. The cores of these rotational vortices are 
normal, or color- flavor unlocked [2l| . If one were to con- 
sider the magnetic field in these cores, and incorporate 
boundary conditions between the cores and the bulk mat- 
ter, one would find a difference between field strengths 
inside and outside of the cores. This difference would 
be responsible for creating a sufficient repulsive force be- 
tween vortices allowing them to drag the magnetic field 
as they move outwards. 

Although 22] found this difference to be small, it is 
sensitive to the value used for the QCD coupling con- 
stant. For example (cf. (22}; ecm - 3.4), with QCD cou- 
pling constants in the range of g 2 /4ir ~ 0.1 to 1, the 
amount of field expelled is 1% to 0.1%, for the case of 
an abrupt transition regionj49j]. As we are considering 
magnetar-strength magnetic fields on the order of up to 
~ 10 16 G at the star's surface, it would then seem reason- 
able that a significant amount in the center is expelled 
from the CFL matter into the rotational vortices. A more 
detailed study, including the running of the coupling con- 
stants [23[ , shows that significant inter- vortex forces can 
exist, but again are sensitive to the masses of the gauge 
fields. At this point the value of the QCD coupling con- 



stant is not well known, and so it is difficult to determine 
whether the inter-vortex force is sufficiently strong. In 
this paper we hypothesize that it is, and note how well 
our model matches observations, and leave further exam- 
ination of the topic for future work. We also note that 
the (possibly more applicable) spin-1 color superconduct- 
ing phases [2J] that do completely exhibit the Meissner 
effect [Hi. 

The total number of rotational vortices at any given 
time is N v = kQ, with k being the vortex circulation and 
f2 the star's spin angular frequency. Hence, as such a star 
spins-down the number of vortices decrease. They do so 
by being forced radially outward [25[ and upon reach- 
ing the surface are expelled from the star. The magnetic 
field, which is pinned to the vortices, is then also expelled 
and the subsequent mag netic reconnection leads to the 
production of X-rays [261 ] - Since the star's spin-down rate 
is proportional to the magnetic field strength, when vor- 
tices are expelled, along with the magnetic flux contained 
by them, the spin-down rate also decreases. Thus, the 
spin-down rate and magnetic field strength of a CFLQS 
become entirely coupled. While, this type of model has 
been proposed in neutron stars [27], in that setting, in- 
teractions between proton and neutron vortices prevent 
the clean coupling between magnetic field strength and 
rotation period. In the CFLQS setting there is only one 
type of vortex. Moreover, the presence of the crust on 
neutron stars further inhibits any clean expulsion of the 
interior magnetic field. In contrast, the vortex expulsion 
process from a CFLQS is very efficient at releasing X-ray 
photons from the stellar surface, as it possesses little or 
no crustal material. As such, the magnetic field is readily 
expelled from the star's interior where it is then able to 
decay through reconnection. 

The X-ray luminosity from this vortex expulsion pro- 
cess is given by [28[ 

L x ~ 2.01 x 10 34 erg s" 1 r) X ,o.i Pl n , (1) 

where ?yx,o.i is the reconnection efficiency parameter in 
units of 0.1 and P the spin-down rate in units of 10~ u 
s s _1 . An estimate for the latter as well as the magnetic 
field evolution, both derived from vortex expulsion, is 

where Pq denotes the rotational period of the star at 
birth, P>o the magnetic field at birth, and r is a relaxation 
time given by 

(3) 

One can see from equations (jT| &; [2j that the X-ray emis- 
sion decreases as the star spins down. This X-ray emis- 
sion is produced on the surface of the CFLQS, which 
alters its thermal evolution. 
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III. THERMAL EVOLUTION 

To study this numerically, the general relativistic equa- 
tions of energy balance and thermal energy transport 
need to be solved. These equations are given by (G — 
c=l) 



dm 

d(Te^) 
dm 



P\J\ — 2m/r 



24> 



d(Te^) 
'~dt~ 



16tt 2 



■(4) 



(5) 



l npy/l — 2m/r ' 

respectively @. Here, r is the distance from the center 
of the star, m(r) is the mass, p{r) is the energy density, 
T(r,t) is the temperature, l(r,t) is the luminosity, </>(r) is 
the gravitational potential, e v (r, T) is the neutrino emis- 
sivity, c v (r, T) is the specific heat, and n(r, T) is the ther- 
mal conductivity Q . The boundary conditions of @ and 
([5]) are determined by the luminosity at the stellar center 
and at the surface. The luminosity vanishes at the stellar 
center since there is no heat flux there. At the surface, 
the luminosity is defined by the relationship between the 
mantle temperature and the temperature outside of the 
star HI]. 

Heating/cooling mechanisms used in our work, other 
than those included in the minimal cooling scenario, 
are the quark direct Urea process, emissivity of which 
is on the order of 10 26 erg s" 1 cm -3 , and the quark 
modified Urea and Bremsstrahlung, which are of order 
10 19 erg s _1 cm -3 . Due to quark pairing in the CFL 
state, however, the direct Urea process is suppressed by a 
factor e _A / T and the modified Urea and Bremsstrahlung 
by a factor e~ 2A / T for T <T C , where A is the gap param- 
eter for the CFL phase and T c is the critical temperature 
below which strange matter undergoes a phase transi- 
tion into CFL matter. In the case of color-flavor and 
color-spin locking it was shown that the critical tempera- 
ture, where the condensate melts, deviates from the BCS 
behavior. In the CFL case, the transition temperature 
is a factor 2 1 / 3 larger than the one would expect from 
BCS theory ([30j]). Here, we assume the validity of the 
BCS relation for the gap A = A x yjl - (T/T c ) 2 , with 
the critical temperature given by T c ~ 0.57A ; A is 
the magnitude of the zero-temperature gap at the Fermi 
surface. We point out that, because vortex expulsion 
becomes dominant at relatively early times, when it is 
included in the calculations the exact value of the crit- 
ical temperature has a negligle effect on the long-term 
temperature evolution. 

In this work the massless Goldstone bosons, due to 
the breaking of baryon number, were not included. In the 
work by 3l| it was shown that at later stages in the ther- 
mal evolution these may become dominant. While [32[ 
also confirm this, both studies look at only neutrino emis- 
sion channels. Further studies including photon and neu- 
trino emission channels [U H3| conclude that the cooling 
time of CFL stars is similar to that of ordinary neutron 
stars. As such, we expect that the presence of massless 



Goldstone bosons would not significantly change the size 
of the shaded band in figure [TJ nor our conclusions. 

The heat produced by vortex expulsion occurs in an 
emission region just above the star's surface. In this re- 
gion, the energy released by the magnetic field decaying 
after it has been expelled from the star's interior, is de- 
posited. An emissivity for this process can be calculated 
from the energy per unit time (cf. Eq. [TJ per unit vol- 
ume, where the volume of interest is a shell surround- 
ing the star. The width of the shell is estimated to be 
the minimum length of a vortex still inside the star, just 
before it is finally expelled. Such a vortex would be a 
distance of approximately the inter-vortex spacing away 
from the surface of the star. The shell width is then 



AR = 8.42xl0" 3 km 



P\ 1/4 / R \ 

Is) v l0kn v 



1/2 



/ 300MeV \ 
{ fx/3 ) 



1/4 



( 6 ) 

where \i is the average chemical potential throughout the 
star. In equation (J6j) it can be seen that the thickness 
of the heating layer depends weakly on the star's spin- 
period and density. This implies that the heating layer 
due to vortex expulsion will not vary significantly from 
one star to another. 



IV. RESULTS 

We have solved equations ^ and (J5J) numerically for 
models whose structure (composition and bulk prop- 
erties) were computed from the Tolman-Oppenheimer- 
Volkoff equations [l|, HJ . The equation of state used for 
CFL matter was the MIT bag model with massive strange 
quarks (cf. for example eqn. 20 in 35]). The parame- 
ters used for the equation of state are m s = 150 MeV, 
B 1 / 4 = 145 MeV, and A as described above. The re- 
sults are shown in Fig. [TJ where the redshifted surface 
temperature evolutions for various types of cooling sce- 
narios are plotted. These scenarios are CFLQSs (CFL 
stars with superconductivity and the resulting vortex ex- 
pulsion), and in the shaded region in Fig. [TJ CFL stars 
(with vortex expulsion intentionally left out), uds stars 
(strange quark stars without pairing), and neutron stars 
2]. The CFLQS birth spin-periods and magnetic field 
strengths were chosen such that the currently observed 
SGR/AXP/XDIN values of spin-period, spin-down rate, 
magnetic field strength, and luminosity would be con- 
sistent with the values derived from vortex expulsion (cf. 
Eq. [5]) . In other words by constraining two parameters in 
our model with observations, the spin-period, spin-down 
rate, magnetic field strength, and luminosity all become 
self-consistent. The observed data is taken from Table [TJ 
where the ages and temperatures of SGRs, AXPs, and 
XDINs are listed. 

Figure [TJ shows that CFL quark stars without vortex 
expulsion cool down too rapidly to agree with observed 
SGR/AXP/XDIN data. One might expect that the sup- 
pression processes that contribute to cooling, due to CFL 
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FIG. 1: Redshifted surface temperature evolution for CFLQSs 
(CFL stars with superconductivity and the resulting vor- 
tex expulsion). The initial spin-periods and magnetic field 
strengths are parameters and were chosen such that (in our 
model) the evolution of spin-period, spin-down rate, magnetic 
field strength, and luminosity would be consistent with obser- 
vations. The shaded region indicates calculations of cooling 
scenarios for various types of neutron stars, uds stars (strange 
quark stars without pairing), and CFL stars with vortex ex- 
pulsion intentionally left out. All calculations were done with 
stellar masses of 1.4 Mq and radii of 10.5 km. The observed 
data is listed in Table [I] 

pairing, would keep these stars hotter for longer, but pair- 
ing also changes the specific heat capacity, resulting in 
enhanced cooling. Standard neutron star cooling pro- 
cesses also lead to stars with temperatures much lower 
(albeit warmer than uds and CFL stars without vortex 
expulsion) than values observed for SGRs/AXPs/XDINs. 

Emission due to vortex expulsion dominates all other 
processes except during the first few minutes. However, 
we also considered the possibility that a neutron star un- 
dergoes a phase transition to a CFLQS after a delay. 
Observational motivations for considering a delay are; (i) 
superluminous supernovae and hypernovae (36|; (ii) the 
large discrepancy between SGR/AXP ages are their pro- 
genitor supernova remnants (37j . Physically, a delay may 
be necessary if one considers the following; i) the time 
needed for a newly born compact object to spin-down 
sufficiently such that the center reaches nuclear densities 
[38j | , ii) the time for the temperature to reach the super- 
conducting critical value (T c ) and the resulting vortex 
lattice to form, iii) the nucleation time for strangelets to 
start fusing together (39j . 

Using the quark-nova model (4(j the star is reheated 
to roughly 10 11 K following the transition to CFL mat- 
ter. This energy is from the release of gravitational en- 
ergy during the collapse as well as the latent heat re- 
leased when converting from hadronic to stable strange 
quark matter [§3, EO, EH . Within the framework of the 



FIG. 2: Thermal evolution of a neutron star undergoing a 
phase transition into a CFLQS after a delay. This delay is 
estimated to be the time needed for the central density of a 
neutron star to reach the critical value at which CFL matter 
is favored (~ 100 years). Values of the birth temperature for 
the resultant CFLQS are of order 10 11 A, and were estimated 
using the quark-nova model [53]. The CFLQS birth periods 
and magnetic fields are as indicated. 



quark-nova scenario the neutron stars with typical values 
for spin-period and magnetic field, and masses greater 
than 1.5 solar masses, are most likely to make the tran- 
sition to quark deconfinement and the subsequent color- 
superconducting phase. Upon this transition the mag- 
netic field may be amplified via color ferromagnetism to 
roughly 10 15 G [42j . although there remains the possibil- 
ity of strong fields from flux freezing during gravitational 
collapse and/or a dynamo at stellar birth. We choose 
initial values for the CFLQS 's magnetic field strength 
accordingly. The results are shown in Fig.[2]for CFLQSs 
born after a delay of tqn = 100 years. 



V. CONCLUSIONS 

We have computed the thermal evolution for various 
types of strange quark stars and found that CFL quark 
stars, possessing a rotational-vortex lattice, are in good 
agreement with observed temperatures of SGRs, AXPs, 
and XDINs. Our model is applicable to any star made 
of a three-flavor color superconducting phase that ex- 
hibits at least a partial Meissner effect (eg. spin-0 [22] or 
spin-1 [l9|). In our model, the CFL star spins-down as 
a result of magnetic braking and expels vortices, which 
contain magnetic flux, thus decreasing the magnetic field 
strength in the magnetosphere, resulting in a lower spin- 
down rate. Hence, the magnetic field strength and spin 
period become coupled. By including emission from 
vortex expulsion in our relativistic cooling calculations, 
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we have found that the unusually high temperatures of 
SGRs, AXPs, and XDINs can be predicted. 

Other distinct signitures of CFL matter in stars in- 
clude a photon fireball [43[ , which is of importance to ex- 
plosive astophysics. However, it is only relevant during 
the earliest (ie. birth) stages of a CFLQS. For long-term 
cooling our findings indicate that vortex expulsion is the 
dominate emission mechanism. 

From previous papers [HI, [28| we have also shown that 
by using our model for a CFL quark star the evolu- 
tion of the spin-period, spin-down rate, and magnetic 
field strength can also be predicted for SGRs, AXPs, 
and XDINs. The long-term evolution of these proper- 
ties suggests an ancestral linkage between SGRs, AXPs, 
and XDINs. We also note a paper which confirms 
that the birth statistics of SGRs and AXPs are consistent 
with the number of observed XDINs. Finally, this study 
suggests that a delay time between supernova and quark- 
nova events (tqn) is possible. If the estimated delay of 
100 years is correct, then the density at which quark mat- 
ter deconfines can be calculated to be roughly five times 



the nuclear saturation density [38j . The quark- nova delay 
can be measured by considering the discrepancy between 
SGR/AXP ages are their progenitor supernova remnants. 
Provided one had accurate measurements of this age dif- 
ference, a more precise value for the density at which 
quark matter deconfines could be inferred. 
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